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Abstract 

We explore the possibility that the observed baryon asymmetry of the universe 
is the result of an earlier phase transition in which an extended gauge sector breaks 
down into the SU(3)c x SU(2)z x U(1)y of the Standard Model. Our proto- 
typical example is the Topflavor model, in which there is a separate SU(2); for 
the third generation from the SU(2)2 felt by the first two generations. We show 
that the breakdown of SU(2); x SU(2)2 — SU (2)z results in lepton number being 
asymmetrically distributed through-out the three families, and provided the SM 
electroweak phase transition is not strongly first order, results in a non-zero baryon 
number, which for parameter choices that can be explored at the LHC, may explain 
the observed baryon asymmetry. 

*jshu@theory.uchicago.edu 
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1 Introduction 


The origin of the observed baryon asymmetry of the universe is one of the most fundamen- 
tal open questions in particle physics and cosmology. Numerically, the baryon-to-entropy 
ratio is ng 

n=- = 9.2406 x 1074, (1) 


where ng = n, — nz and s are the baryon number density and entropy, respectively. 
Famously, Sakharov has shown [Í| that in order to generate a baryon asymmetry in an 
initially baryon-symmetric universe, there must be: (1) baryon number (B) violation; (2) 
C and CP violation; and (3) a departure from thermal equilibrium. These requirements, 
especially in light of null results obtained by experimental searches for B-violation (for 
example, proton decay), are difficult to realize in the framework of an electroweak (EW) 
scale model, prompting attention to high scale mechanisms such as GUT baryogenesis 
[4 [5], leptogenesis [6] [7] and the Affleck-Dine scenario [8]. Such mechanisms can be 
viable, but are difficult to test experimentally. 

However, the Standard Model (SM) already contains a way to reconcile large baryon- 
number violation at the EW scale with the lack of experimental evidence for such inter- 
actions at low scales. As noted by Kuzmin, Rubakov and Shaposhnikov [9], the baryon- 
number violation induced by sphaleron transitions between two different N-vacua of the 
SU (2), electroweak interaction are unsuppressed at temperatures larger than the EW 
scale, but are so highly suppressed at low energies as to be essentially irrelevant. This 
idea of electroweak baryogenesis (EWBG) is very attractive, and could have been re- 
alized within the SM. EWBG makes very specific demands of the electroweak symmetry- 
breaking (EWSB) sector of the theory, and leads to observable consequences at future 
colliders. For example, the need for B-violating interactions to be out of equilibrium at 
the time of the phase transition requires that the electroweak phase transition is strongly 
first order, with 


(2) 


This in turn puts constraints on the Higgs potential, and demands a small Higgs quartic 
(and hence Higgs mass). 

Electroweak baryogenesis in the SM is thus ruled out; the LEP-II bound on the 
SM Higgs mass [I], m?™” > 115 GeV is incompatible with a the requirement that the 
sphaleron processes be out of equilibrium at the phase transition, Eq (2). While this ren- 
ders the SM itself unable to produce the asymmetry, it illustrates the fact that experiments 
at the TeV scale are able to directly probe EWBG. Extensions of the SM which include 
new physics at the EW scald] can in fact re-open the possibility of EW baryogenesis by 
introducing additional bosonic fields [2J 03, [14] |15} [16] as in the Minimal Supersymmetric 


‘In addition to new fields, a greatly increased Hubble expansion{27| can re-open some of 
the Higgs mass parameter space. Such models usually require the very fast Hubble expansion to only 


happen close to the electroweak scale. 


Standard Model (MSSM) or fermions strongly coupled to 
the Higgs [32]. Eventually, future colliders such as the Large Hardron Collider (LHC) 
and possibly International Linear Collider (ILC) will unravel the nature of EWSB, and 
should provide a better understanding as to the nature of the electroweak phase transi- 
tion (EWPT). If it proves first order, this will be a crucial piece of indirect evidence for 
EWBG. If it does not, it will raise the interesting question: Does the idea of electroweak 
baryogenesis have to be abandoned? 

While the lack of a first order EWPT would strongly disfavor electroweak baryogenesis, 
it may be that the basic picture of baryon-number violating processes arising from non- 
perturbative gauge theory dynamics, impotent at low energies but unsuppressed at the 
TeV scale, can survive. Since many extensions of the SM predict the existence of new 
non-Abelian gauge symmetries at the TeV scale, it may be that some theories already 
addressing unrelated problems may in fact contain the ingredients necessary for “electro- 
weak style” baryogenesis. A new gauge group G could break down somewhere above the 
EWPT, and generate an asymmetry through its own strongly first order phase transition. 

This is a novel idea, but one which is not entirely straight-forward to realize in a 
realistic setting. Below this new phase transition, the theory is still SU(2); x U(1)y 
symmetric, and the usual EW sphalerons are unsuppressed. They will try to wash-out 
any generated baryon asymmetry unless B — L # 0. In principle, one could arrange the 
representations of the SM matter under the new gauge symmetry G such that B — L is 
not preserved due to the nonzero mixed anomaly G-G-U(1)s—1, but this kind of non- 
universal assignment of representations is also somewhat delicate. The requirement that 
G-G-U(1)y gauge anomalies cancel will in general require more chiral fermions charged 
under G and the SM. Thus, we will not consider this possibility in detail. Instead, we 
employ a more subtle mechanism in which G acts only on the third generation fermions, 
and thus B — L is conserved (within a generation). The breakdown of G can thus be 
accompanied by production of an asymmetry among the third family quarks and leptons. 
Baryon number will quickly equilibriate among the three families because of the large 
quark masses and mixings, but each lepton family number will remain distinct because 
of the tiny neutrino masses. Above the EWPT, the EW sphalerons will result in B = 0 
and Lı = Ly = —L3/2 # 0. Provided the EW sphalerons are in equilibrium during the 
EWPT, as the fermion masses turn on from EW symmetry-breaking, a non-zero (though 
diluted) baryon number will result [83} [34]. So in fact this mechanism requires that the 
EWPT not be strongly first order. 

The outline of this paper is as follows. In Section P| we review how a baryon asymme- 
try can be generated even when B—L = 0, provided there is a non-zero asymmetry in the 
third family lepton number, and how this eventually translates into a baryon asymmetry 
after the EWPT. We apply this mechanism to the “Topflavor” model which is known 
to contain non-perturbative interactions which violate baryon- and lepton-number in the 
third family [86]. In fact, we find that it is possible to generate a baryon asymmetry of the 
right magnitude. The Topflavor model, phase transition, CP violating sources, diffusion 
equations and calculated baryon number density are considered in Section [B] Our results 
show that the right baryon number density can be generated for parameters that would 
render this model testable at the LHC. In Section [4] we present our conclusions. 


2 B froma Family Asymmetric Distribution of L 


In this section we review the mechanism by which an initial condition that has B = L = 0 
can nonetheless result in a non-zero baryon number density provided Ly = Ly = —L3/24 
0 [83} [34]. We will show how the specific example of the Topflavor model can generate 
these initial conditions through the non-perturbative dynamics of its phase transition (in 
a very similar way to that in which baryon number is generated in a traditional EWBG 
scenario) in Section Thus, for now we assume that the initial phase transition has 
generated La = B = A and Lı = Lə = 0. The unsuppressed EW sphalerons will rapidly 
evolve this into a state with B = L = 0 but Lı = Lə = —A/3 and L3 = +2A/3. We 
now study what happens as the universe moves through the EWPT, and show that the 
non-zero T lepton mass will result in B ~ 10~-°A. Our discussion closely follows that of 
Ref. [34]. 

The smallness of observed neutrino masses indicates that lepton-number violation is 
out of equilibrium at the EW scale. Thus, each lepton flavor has a seperate chemical 
potential u; with ¿ = 1,2,3. This is in contrast to baryons, because the large quark 
masses and mixings keep baryon flavor violation in equilibrium, and thus the quarks are 
described by a single chemical potential u. We consider the SM matter consisting of three 
families each of which consists of two quarks (an up-type and down-type) with masses Mg, 
a charged lepton of mass m, and a neutrino which for our purposes can be approximated 
as massless. The free energy per unit volume for the system in equilibrium at temperature 


T is given by 
6 3 


F= 6X F (ma, 4) + XOF (mu, pi) + F(0, 13)] (3) 


i=1 i=1 
where the SU(2); gauge interactions maintain equilibrium between the charged lepton 


and its neutrino and the up- and down-type quarks of a given generation. The free energy 
density for a (single helicity of a) fermion of mass m and chemical potential u is given by 


3 


F(m, u) = -T / Saal + e7 HO 4 In(1 Ag ET], (4) 


where K; is the spatial momentum of the fermion, and E = \/K? + m? is its energy. At 
high temperatures, T >> m, u, this may be approximated as, 


1 3 m?’ 
F (m, u) ~ F (m, 0) — peT (1 -= t) . (5) 


The (individual) leptonic and baryonic number densities may written 


d 1 
Li = iy Hi D| © —5uT’ bi 
gu eet ED) HiT Pi, (6) 
de 1 
— N 2 


i=1 


where 
6 2 


z 3 Mg; = 1 mi, 
i=1 


Electroweak Sphaleron transitions violate X` L; and B, but preserve the three combi- 
nations A; = L; — B/3. In terms of the chemical potentials these are 
1 uT? uT? 


A; = L; — =B x —a- 
I= g 9° 2 


Gi, (9) 


We can invert the above relations to obtain each u; in terms of the quark chemical potential 
L, temperature T, and the conserved value of the corresponding A,;. Effectively, the EW 
sphalerons convert nine quarks and one lepton of each family into nothing. In thermal 
equilibrium, this leads to the relation u = — >>, 4;/9. Using this fact, together with the 
three conservation equations Eq (9), allows us to express the quark chemical potential in 
terms of the values of the A,, 


(REDE a 


i=1 


which can be combined with Eq (6) to obtain the final baryon number density 


4 
—(B-L B-L#0 


4 N i 
-Br Din Aire 


The first of these results is the familiar relationship applicable to theories that directly 
generate a non-zero B — L (such as leptogenesis) and indicates that in such theories 
primordial B cannot be completely washed out, and a primordial L will be converted into 
B by EW sphalerons. The second result shows how in a theory with B = L = 0 but 
the individual A; non-zero, the turn on of the charged lepton masses will also generate a 
non-zero B. In the scenario we are considering, with initially B = 0 and L3 = 2A/3, and 
taking the freeze-out temperature to be the close to the EW scale, the resulting baryon 
number is diluted to about B ~ 10~°A P, 33) BA. 

Since the dilution factor plays a relevant role in our work, let us expand on its origin: 
To compute the above quoted dilution factor, we have assumed a second order phase 
electroweak phase transition. Under this condition, the sphaleron processes will remain 
in equilibrium until the weak spahleron rate is of the order of the expansion rate of the 
Universe. The departure from equilibrium therefore occurs at the freeze-out temperature 
Tr, such that v(Tp)/Tp ~ 1. Using the relation m,(T) ~ h-/V2v(T), and the condition 
u(Tr)/Tr = 1, we get that the final baryon number is approximately given by 


(11) 
B-L=0. 


2 
4 Al mw 16x 10-8A. (12) 
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3 Lepton Number Generation in the Topflavor Model 


3.1 The Topflavor Model 


The gauge extension of the SM that we consider is based on the gauge group SU(3), x 
SU(2), x SU(2)2 x U(1)y. While the SU(3),. and U(1)y subgroups remain the same as 
those of the SM, the SU(2), group of SM is expanded to a larger SU(2); x SU(2)2 ina 
flavor dependent way. The fermion content of the model is identical to SM. Under the new 
SU(2), and SU(2)2 groups, the doublets of the third generation transform as doublets 
under SU(2), and singlets under SU(2)2, while the first and second generation doublets 
are singlets of SU(2), and doublets of SU(2)z. Thus, their SU(2); x SU(2)2 x U(1)y 


quantum numbers are 


Q? = (2, 1)1/6, qr = (1, 2) 1/6, 
D = (2,1), D? = (1, 2)-1/2- (13) 


After symmetry breaking, the Standard Model SU(2); group emerges as the unbroken 
diagonal subgroup of SU(2); x SU(2)z. The corresponding SU(2); gauge coupling is 


N92 


g, = ——_.. 
Vg +g 


This relation implies that when one of the gauge couplings becomes large, the other one 
approaches gz from above, and thus both gı and g2 are necessarily larger than gz. Thus, 
a convenient parameterization is given by, 


(14) 


gL 
= 22 15 
g2 = COS , ( ) 


— IL 
gı = sino’ 
in terms of an angle ¢. We will use a simplified notation s = sin@ and c = cos ¢ below. 

The symmetry breaking of the extended gauge group, SU(2), x SU(2)2 to the SM weak 
gauge group SU(2),; is accomplished by introducing a vacuum expectation value (VEV) 
to a scalar field ©, which transforms as a bidoublet (2,2) under the extended gauge group 
transformations. After the SU(2), x SU(2)2 breaking, © can be decomposed under the 
residual diagonal SU(2); symmetry into a complex singlet o and a complex triplet 7 (half 
of which is eaten by the Z’ and W’s), 


o + T3 Jot, 


1 
y= 
2 J2r_ 0-7 


(16) 


We introduce the scalar potential, 


1 ~ 
Va = mB? + Al(EE)P + NIE + (-3DEE) + XEDE? + hc), (17) 


where D and À are complex parameters and we use a notation that supresses the gauge 
indices: (UX) = Yyp™-g€acépg and |X|? = Tr(HtH). For appropriate choices of parameters, 
this potential results in the VEV, 
Ly 
(Za) = Zuo eld (18) 


which will generally be complex, and will be solution of the equations, 


D 2100 D* —2i09 — 2 
E rA (19) 
N+ A + Ae?* + \*e7 260 
1 ~2D* + A'u? 
ło = —-acos Re ea (20) 
4 —2D + dug 


Choosing some representative parameters, taking m = 200 GeV, D = 5 x 108 et GeV?, 
and à = A = Aà = 0.05, we obtain a zero temperature VEV described by uo œ 2.7 
TeV and 6) ~ —0.7. This particular set of parameters has been chosen with small quartic 
interactions in order to have a first order phase transition, with the dimensionful quantities 
arranged such that the SU(2) x SU(2) breaking scale is of order TeV. Precision electroweak 
constraints have been extensively considered in the literature [35], and typically require 
Ug > a few TeV. Requiring that the extended instantons of the strongly coupled SU(2) 
do not mediate unacceptably large proton decay further requires the gauge couplings 
to satisfy sin? @ > 0.2. We will illustrate our results with the representative point chosen 
above, and s? = 0.4, in order to be consistent with all constraints. 

The fermion doublets of either SU (2); or SU (2)z transform as doublets under SU(2),. 
In addition, there is a SU(2),, triplet of heavy gauge bosons from the breaking. We denote 
the neutral and charged heavy gauge bosons as Z’ and W’~. Their masses are degenerate 
and given by 


Mwy = Mz = (9? + g3)u? = Er. (21) 

At the electroweak scale, v ~ 174 GeV, the remaining SU (2); x U(1)y electroweak 
symmetry is broken to U(1)em as in the SM. This is accomplished by giving a VEV to one 
or more Higgs boson doublets. There are two possible representations for the Higgs boson 
under SU(2); or SU (2)z, either (2,1)41/2 or (1,2)41/2. We will focus on the first case 
(sometimes called the heavy case) as it motivates the large third family fermion masses as 
they are the only family whose Yukawa interactions are SU(2); x SU(2)2 gauge-invariant. 
In a non-supersymmetric theory, a single Y = +1/2 Higgs doublet H, suffices, but our 
results are largely unchanged if we generalize to a Y = +1/2 pair of doublets, H,, Ha 
instead. In order to connect more easily with the supersymmetric case, we consider the 
case with two Higgs doublets below. 

The Yukawa couplings for the first two generations can be generated by adding an 
additional “spectator” Higgs-like doublet H/ (in a supersymmetric theory it would be 
a pair of doublets including H/,) charged under SU(2)2. They couple to the first two 
generations via Yukawa couplings and mix (slightly) with the regular Higgs(es) via inter- 
actions such as A, HH}. The small Yukawa couplings for the first two generations can 


7 


be naturally obtained from the small mixing angle between H’ and H. The full set of this 
kind of gauge-invariant =-H-H’ interactions between Higgs(es), spectator Higgs(es), and 
X include, 


A\H, SH +h.c;  AgHDH' + h.c, 

A’ H,.S' Hi, + h.c; A, H{™' H! + h.c, 

cou* H! DH! + h.c; ctu! H! E HH + hc, 
cy HSH) + h.c; chu! HST HT + hc (22) 


These interactions will be important below, because they (indirectly) drive the Topflavor 
phase transition’s production of Ls. 


3.2 The Phase Transition 


The details of the SU(2), x SU(2)2 — SU(2)z phase transition depend on the finite 
temperature effective potential. In order to generate a baryon asymmetry in the Topflavor 
model, the processes which violate baryon- and lepton-number must be out of equilibrium 
at the time of the phase transition. We are interested in the regime where gı > go, for 
which we can approximate the sphalerons associated with the SU (2); x SU(2)2 symmetry 
breaking as purely arising from SU(2),. Their rate, 


7 
u 
Eph =~ &K T6 Exp [—Espn/T] , (23) 


(with Espn = 4Tu/gı and « a dimensionless parameter of order one [87]) must be much 
less than the Hubble expansion, H ~ gT? mp. Assuming «s ~ 1, gą ~ 100 and 
ue Z Te ~ TeV, this requires, 
Uc 
nt 
and for s? ~ 0.4, we should havd tif Ta Z 2.5. 
The potential for X, Eq. @7), can be expanded to, 


Z 2.5, (24) 


Vo(o) = 


A+ 1 A 
=m ||? + ——|o|* + | —=Do? + —o"|o"| + h.c. 2 
5 |o| A |o| í ma |o*| c (25) 


where for brevity we have not shown the terms involving the triplet component 7. At one 
loop, there are both temperature-dependent and temperature-independent corrections to 
the potential, 


V(u,9,T) = Vo(u,0,0)+ Vi(u, 0,0) +Vi(u, 8, T) . (26) 
2Note that such a strong first order phase transition may provide a strong signature in gravitational 


waves detectable at the planned space interferometer, LISA. We will pursue this idea 
in a separate paper|45). 


We consider the limit where the gauge couplings of the SU(2)’s are much larger than 
that of the © self-interactions À and thus approximate the complete one-loop corrections 
by those from the gauge bosons, Z’ and W’. We renormalize parameters according to a 
scheme that absorbs the correction to the position of the zero-temperature minimum, and 
find, 


2 
nh Oy = 2 2 ela igs 2 ee (27) 
WBD T Gage gaz J | L008 ue 2 “o 


and finite temperatue correction, 


Vi(u,0,T) = ies petal (28) 


with 


io) = i. dx- x {log (2 —e Vv een Pear! — log(1 — en (29) 


T 0 i 


where g; = 9. The Debye screening effect on the longitudinal modes of massive gauge 
bosons is neglected because their corrections gT is small compared to their masses in- 
duced by the Higgs vev|47]. Note that the one loop corrections from the gauge sector 
depend only on the magnitude of the VEV u and not on its phase. 

For the sample parameters chosen above, m = 200 GeV, D = 5 x 10° e GeV?, 
à = X =X=0.05,and s? = 0.4, we find that the critical temperature for these parameters 
is 7, ~ 840 GeV, and the VEV at T, is described by ue ~ 2.7 TeV, 0e ~ —0.7, indicating a 
first order phase transition that is easily strong enough. In Figure [I we plot the effective 
potential for several choices of temperature. 


3.2.1 Bubble Profile and Evolution 


At T., the rate for nucleation of bubbles with (7) 4 0 becomes large, and the bubbles 
expand to fill the universe with the true vaccuum. In this subsection we make some rough 
estimates of the properties of the nucleated bubbles, which are pertinent to the eventual 
generation of baryon asymmetry as they provide the out-of-equilibrium dynamics which 
results in lepton-number being unequally distributed through-out the three generations. 
We will simplify the treatment by considering the phase transition as a quasi-equilibrium 
process such that the temperature is varying slowly enough that the properties of the 
nucleating bubbles can be obtained by studying fixed temperature solutions at T ~ Te. 
Assuming that the nucleated bubbles have spherical symmetry [88], the Euclidean 


action of the configuration becomed}| 


S(T) = 4r / a r?f (Vo) +V((0),7.)} (30) 


3We use the O(3) approximation since the nucleation temperature Th ~ Te > 1/(2R) [£6]. 


Figure 1: The effective potential as a function of the magnitude of the VEV for three 
different choices of temperature. The phase of the VEV at each point is chosen as the 


solution of the Equations of motion for that value of the magnitude. 


where (ø) is a function of radius r that describes the configuration. The transitions will be 
predominantly mediated by the solutions which minimize this action, give by the solutions 
of the equations, 


u 22du ôV 


dr? ` rdr Ou 
dO 2d0 ôV 
. E r | 60 ey 
subject to the boundary conditions, 
a = (} ul, œ = 0 
dr | o 
d 
a 0; Olo = u= - (32) 
dr |o 


These coupled differential equations are somewhat difficult to solve. Instead of looking 
for detailed solutions, we will use an ansatz for the profile, and a variational approach 


10 


to determine the parameters that describe it. We write the solutions in the form of two 
“kinks” with the proper asymmptotic behavior, 


= 
pm 
3 
wy 
| 


= [1 — Tanh (a (r — R))] 


Oe — bu- 
olr) = E 


[1 — Tanh (a (r — R))] (33) 
where the radius of the bubble is ~ R and the width of the bubble wall is ~ 1/a. In 
addition to imposing the form of the solution, we also assume that the bubble width is 
the same for u and for 0 (which we expect to be approximately true). We determine a 
numerically by plugging the solutions of Eq. (83) into the action Eq. (80) and finding 
the value of œ which minimizes the action. For our standard parameter choice we find 
a ~ T/2 and Ly = 1/a ~ 2/T. The profiles are plotted in Fig 2] We expect that for 
large bubbles, the details will become independent of R, which in fact proves to be true. 


3 -0.45 
2.5 -0.5 
2 -0.55 
3 D 
1.5 -0.6 
1 -0.65 
0.5 -0.7 
o 2 4 6 8 10 12 14 i6 18 07” 
r (TeV) 


Figure 2: The bubble profile at the critical temperature T, for R = 10 TeV’. 


Under our quasi-equilibrium assumption, the expansion of the bubble is driven by the 
fact that the gauge bosons acquire masses inside the bubble, and thus the free energy is 
minimized for large bubbles [47] [39]. By analogy with the SM EW phase transition, we 
estimate the bubble wall velocity vw ~ 0.05, though we find that our final results are very 
insensitive to the precise value of vy. 


11 


3.3 Diffusion Equations in the Topflavor Model 


We now compute the prediction for the L3 generated during the transition in which 
the Topflavor model breaks down to the SM. The underlying picture is similar to the 
standard EWBG picture in the SM (or MSSM). The bubble of true vacuum is expanding, 
and generates chiral charge through the C’P-violating interaction of the plasma with the 
bubble wall. In the specific case of Topflavor, the particles which interact strongly with 
the wall are the Higgses, through the interactions in Eq (22). These charges diffuse freely 
in the unbroken phase and are converted into B and L3 by a combination of the Yukawa 
interactions and the unsuppressed sphalerons. As they pass into the broken phase, they 
are frozen. 

In the limit gı >> g2, we neglect the SU(2)2 sphalerons associated with the first two 
families. The quark Yukawa interactions and the QCD instantons, together with the fact 
that all of the quarks diffuse at approximately the same rate, allows us to constrain the 
light quark densities in terms of the right-handed bottom density b, 


QiL = Qə, =U p= —2D 2 = 299 = —2Cpr = —2b. (34) 


Thus, the species whose densities we will track are the left-handed top and bottom doublet, 
Q = tL +b, the right-handed top t = tr, the right-handed bottom b = bp, the left-handed 
lepton doublet L3 = T; +v,, and the Higgs h = (ht +h? — h3 — h9). We assume that the 
H-H'-X interactions Eq. (22) are fast enough such that the spectator Higgses H’ are kept 
in equilibrium with the Higgs, and thus © = 0, h! = (h/t + h? — h'g — h’°) = h, and we 
do not include the densities © and h’ in the diffusion equations. For relativistic particles 
near equilibrium, we can write the number densities in terms of a chemical potential p 
as ni = kju;T?/6 where k; counts the number of degrees of freedom, 


kg =6; kr = 2; neha 3; kp = 8. (35) 


The diffusion equations will contain the interactions which are fast compared to 
the time scales at which the elements of the plasma are diffusing, 7; = D;/vw, where 
Uy ~ 0.05 is the speed of the bubble wall’s expansion and D; is a diffusion coefficient 
which characterizes the interactions with the background plasma. Typically, one expects 
Dg = D; = D, ~ 6/T and Dz ~ D, œ 110/T [49]. Thus, we consider the processes 
characterized by rate I = Tg. These interactions include the SU (2), sphalerons with rate 
Ty, the QCD instantons with rate [gc¢p, and the top quark Yukawa coupling to the Higgs 
with rate [,. We continue to assume that the sphalerons associated with SU(2)2 can be 
neglected. These rates are estimated to be equal to BOJ, 


27 Ay 
T > TAAS (2) T = 7.4 GeV, (36) 
Teco ~ 16k'aģT = 0.3 GeV, (37) 
Tr, œ 30a,°T = 0.1 GeV, (38) 


where à ~ 1 is the top Yukawa coupling, as(T.) ~ 0.08 is the strong coupling constant, 
and K’ ~ 1 is a dimensionless coefficient. We have evaluated the rates for s? = 0.4 and 
T, = 840 GeV, as is appropriate for our example parameter set. 
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We approximate the bubble as large, and thus treat the problem one-dimensionally, 
with the z-axis perpendicular to the wall, whose location is at z = 0, with the z > 0 side 
in the broken phase. The rates of change of the various densities are described by the 
coupled set of equations4 


Q h t Q t b 
vQ! = DoQ" E p a T J a aG 
vu — DoQ ae kn i OOP ko k T 
Q Ls 
=6f 324 2 
peri ’ 
Q h t Q t b 
=D =: PE a y a 
” Q al ee el ee le u 
Q h t 
Voh — Dih” = -T -2+4 +n, 
h y ko kp, ki Yh 
Q t b 
„b — Dob!” = Sr PERREN: a 
Vwb Q 3 oo 2 ki i ’ 
L 
toby = Dr = -N b2+e], (39) 
ko kr 


where primes denote derivatives with respect to z and yp is the C'P-violating source for 
the Higgs induced by the bubble wall, approximated as a step function, 


Yhu — n (—Ly < z <0) 
Yh = , (40) 
0 (z>0 or z<—-L,) 


whose magnitude we estimate below. 


3.4  C'P-violating Sources from Spontaneous CP violation 


We consider the C’P-violation arising from the spontaneous C'’P-violation associated with 
the phase of the VEV (c). This field couples directly to the EWSB Higgses H, and Ha 
and thus influences their number densities as they scatter off of the bubble wall. As we 
saw in section the phase of (a) varies as one moves from inside the bubble of true 
vacuum to the unbroken phase (see i.e. Figure B), 


Ap + Mu? 


= +s (41) 
—2D + àu? 


1 
Ad = be— ĝÂu=0 = — acs | 


and thus is space-time dependent as the bubble expands. 


“Note that leptons diffuse faster than quarks, and thus the B — L charge density (Q + t + b)/3 — L is 


not zero locally. 
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2 2* 
Q Myg mep O] Miy MwA 


Figure 3: Feynman diagrams for the leading contributions to the self-energies of H„ in 


the background of a spacetime varying (X). 


In computing the value of the C'P-violating source yn, we follow the treatment first 
introduced by Riotto based on the closed time path (CTP) formalism, which 
allows us to capture the main non-equilibrium quantum effects. The CTP formalism 
distinguishes fields with arguments on the positive and negative branches of the closed 
time path. This doubling of fields leads to six different real-time propagators; for a generic 
scalar field œ, 


G3(2,y) = —i(¢(x)¢"(y)), 

G3 (x,y) —i(o'(x)o(y)), 

Gy(z,y) = O(a,y)G3 (x,y) + Oy, 2)GF(a, y), 

Gi(x,y) = 0(x,y)GS(a,y) + O(y, x) GF (x,y), (42) 


which are conveniently written as a matrix G(x, y): 


G(x y) _ G(s, y) -G5 (z, y) , (43) 
G>(z,y) —G*(a,y) 


From the Schwinger-Dyson equations of the path-ordered two-point functions, one obtains 
Xo 
Bii = — f dbs f’ dzo [EZ (X26512, X) - G(X, E$ (e X) 
+G (XE (z,X) -EXE eA- (44) 


The leading contribution to the self energies © y, (x,y) and $p, (x, y) from the inter- 
actions of Eq. are (see Figure B) 


Sp, (£,y) = gl£, y)Ĝn, (x,y) + h(x, y)Gur (2, y), 
Va (z,y) = G(x, y)Gu (x,y) + h(x, y)Guy (x,y), (45) 


where 


gley) = Mal a U) = SUE) + ey Ee) lon Eu) + ciw’="(y)] 
h(x,y) = mig(a)m* aly) = [AE nee Ht (e)] [AE (y) + ATZy)] 

i(z,y) = mMig(x)m7 Gay) = oe a) + cip="(2x)] [e ane ) + ap*X(y)] 
h(c,y) = may (x) wa) = [AE eee )] AE) + ASD(y)] (46) 


The CP-violating part of Eq. (44), evaluated for H„, is the imaginary part of Green 
functions, 


Yh, = -i f ae [- dzo4 [g(z, z) — g(z, x)| Im Gk, (x, z)G§, (z, x) — Gr. (x, 2) Gi, (z x)| 


+ [h(x, z) — h(z, x)| Im GR, (x, z)Gq, (z, x) — Ga, ( x, z) Gq, (z, x) | i (47) 


We can express the scalar Green’s function in terms of the Bose-Einstein distribution and 
the spectral density of a scalar field 


dtk , . 
z = —ik-(a— z 
Gle) = | EE Do melko), (48) 
where the equilibrium distribution functions are: 
IB (w, LL) = NB (w, pis); (50) 


with np(x) = 1/lexp(#/T) — 1]. The spectral density in the limit of small decay width is 
52), 


D = larer a e e y a E 


where w?(|k|) = k? + M?. The mass M and width I are the thermal mass and width, 
respectively. 

The coefficients |g(x, z) — g(z, x)| and [A(x,z)— h(z,x)] can be calculated to the first 
order in the expansion about x = z. 


2i [Je]? — |catt!”] sin [0(e) — 8(2)] u(x)u(2) 
2i [le — leoul?] (£ — 2)" (0,0(a)]w(a) 82) 


When these are inserted in Eq. (47), the fact that the spectral density is isotropic in space, 
implies that only the time component is non-vanishing, and we can make the replacements, 


g(x, z) a g(z, x) 


K 


dedma S nia jae (Fo) wa), e3 
h(a,z)—h(z,2) = 2i [|4]? — |411] (x — z)? (o) u’ (x), (54) 
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and thus, 


i Ab 
Ve = (Fe) u*(2) { [lew ~ leotel”| Zinc, + [l4]? — 141°] Tut} PP) 


lo) k2 f 
aE | dk x {(1 + Refs (wy + Pay)]) om P,P) 
i 0 2m WHWH, i i i į 

+(1 + 2Refng (wa, + iFa, ) Dwn, Tu, waa Tua) 


—2(Im[ng (wa, + iln, )] + Im[ng (wy; + ilu )])G (wa, Tuu, we, Pi) } , (56) 


and the functions J and G are given by 


1 1 a+c 
I(a,b, c,d) = 3 [late +t ay sin E arctan a <] 
1 1 , a—c 
EET sin [2 arctan Ia <] 
(57) 
1 1 a+c 
G(a, b,c, d) = Eo cos E arctan Ta =| 
1 1 a—c 
-aeea ad cos [2 arctan x <| ; 
(58) 


In exactly the same way, we derive, 


7 IN, 
Ve = (He) u (x) { (lew? = |eop'|?] Lin + [|A2l* — |431] Tram} > dea 


The over-all magnitude of the CP-violating source depends sensitively on several param- 
eters which have up until now not played a large role in deriving our ressults. Thus, we 
content ourselves with the order of magnitude estimate based on the sample parameters 
for the ø potential and bubble wall velocity and profile, and assuming the thermal masses 
and widths are roughly TeV, and that the X-H-H’ dimensionful interactiond)| are of order 
TeV. Evaluating the integrals numerically, we find y, ~ 0.01 TeV A@/L,,v,u2, which for 
the choice of sample parameters described above leads to yp ~ 10° GeV“, and acquires 
non-vanishing values only at the position of the bubble wall, where the Higgs fields are 
varying. 


©The thermal widthes of the Higgs and spectator Higgs are dominated by the trilinear interactions 
X-H-H'. 
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The results presented above rely on a method similar to the one used in Ref. 
in the Minimal Supersymmetric Standard Model case. Further improvements to this 
method were performed in Refs. [21], by considering the all-order resummation of the 
Higgs mass insertion efffects. The result was a mild suppression of the results in the case 
of degenerate masses, as well as new relevant terms away from the degenerate mass regime. 
Furthermore, in Ref. [22], based on self-consistency arguments, a more detailed analysis 
of the relation between the CP-violating sources and the currents induced by the Higgs 
fields was performed, leading to the presence of higher order derivatives in the sources, as 
first suggested in Ref. [I7]. The final result of these investigations is a suppression of the 
baryon asymmetry by a factor of a few in the degenerate mass regime compared to the 
one obtained in Ref. [19], as well as a power-law suppression with the masses when they 
move away from the degenerate region. 

One of the weaknesses of the above described work is the lack of a rigorous derivation 
of the sources for the diffusion equation. A derivation of the semiclassical forces in the 
transport equations by means of the dynamics of the two-point function in the Schwinger- 
Keldysh formalism was performed in Refs. [23} [24] 25, [26]. In particular, in Ref. [26] a 
consistent treatment of the fermion mixing effects was performed, and a careful derivation 
of the sources in the diffusion equations was obtained. The final result was smaller 
by a factor of a few to an order of magnitude than the result obtained in Ref. in 
the degenerate mass regime, and a stronger than power-law suppression away from the 
degenerate case. In our work we are interested in an order of magnitude estimate of 
the result for the baryon asymmetry, and for that purpose the results of this section are 
sufficient. However, while an enhancement of the sources by a factor of a few would be 
easy to obtain by a careful choice of the parameters of our model, an enhancement by 
several order of magnitudes woulde prove very difficult. For these reasons, it would be 
interesting to pursue a more general treatment using the techniques of Ref. [26] to make 
a detailed exploration of the full parameter space consistent with the observed baryon 
asymmetry. 


3.5 Results 


We assemble the results and make a prediction for the density of lepton number stored 
in the third family, L3. We have solved the diffusion equations in two ways; the first is 
a brute force numerical solution of the coupled differential equations, while the second 
proceeds by making some simplifying approximations which allow us to determine an 
analytic solution to the diffusion equations. Our analytic solution simplifies the problem 
by assuming that T}, Pacp and (for z < 0) T; are all strong enough that they enforce 
near-equilibrium relations among the number densities of the species participating in the 
interactions. Thus, we have, 


~O, (60) 


and for z < 0, 
j—— —- N _ _ — N 0. 6l 
tne (61) 


Following the treatment of the usual EW case [87], we take linear combinations of the dif- 
fusion equations which are independent of T}, V@cp, and I; and use the equilibrium 
relations to express the remaining two equations in terms of Q and h. The result can be 
expressed as a matrix differential equation, 


we egg hen) (62) 


where we simplfy yp as constant for — Lu < z < 0 and M and N are constant matrices, 
functions of the D’s, k’s, and vy, 


ao [ Pelt) meote) afa F] 
oes o(a) oof) 7 [a 
Ike—ky, __ 9kp+9ki tke 4 _17 
De = Uy ee ee. = Uw 4 = (63) 
Ike Ikn 9 3 
We convert Eq. (62) into a first order differential equation by defining Y Q, hJ”, 


H =W’'. We solve the equations seperately for z < Lẹ (region I), —Ly, < z < 0 (region 
II), and z > 0 (region II) and match W and H across the boundaries. 

In regions I and III, ya = 0, and the solutions are those of the corresponding homo- 
geneous equation, 


Him (2) = Hy (2) Drm = Exp [-M-'Nz] Dım (64) 
Yim(z) = / dx Hy (x) Dim + Cum (65) 


where Cim and Diy are vectors specifying the boundary conditions. The exponential 
terms grow with z, and thus the requirement that U(z — —oo) — 0 requires Cr = 0 and 
the requirement that Ų remain finite as z — +00 requires Dy = 0. 

In region II, the solutions will be the sum of a homogeneous solution with integration 
constants Cy and Dy, and a particular solution, 


ae Ss ae f A an aM : , (66) 
So, 

Hu(2) = Hu (2) Du + HeC) , (67) 

tee). = f de {Hy (2) Pu + Hp (2)} + Cn, (68) 
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and matching this to the solutions in regions I and II determines Dı = Dy, and, 


0 
Dy = f di (2)M3| "| , (69) 
a 0 
—Ly 
Cr = f dxHy (x) Dr , (70) 
0 0 T Yh 
Cm = f dxH yy (x) Diy — [. dx L dyH y (x) Hy (y) Mt ‘ ; (71) 


Note that ¥(z) = Cm for z > 0, so this last expression is in fact the final densities 
produced by the phase transition. 

Since Ly is small, we can derive an approximate form based on the limit Lu —> 0. To 
leading order in Lọ, we have 


0 
Ci = Ci = / dvHy (x) Dr = — La ie an ; (72) 
= 0 


where 


Dy = L,M| |. (73) 


0 


Using the expression for N~', we find the final lepton number density 


12D Yn 
227Vy 


Recalling equation (58) that yn X v,/Ly, we find that our final result is approximately 
independent of the diffusion constants, the bubble wall velocity, and the bubble wall 
width, as long as I’; is fast enough. 
Assembling the results, in Figure [4] we plot the densities normalized to the entropy, 
s ~ 2n?/459,.T?, where gẹ ~ 100. The densities L3, t, and b are determined using the 
relations Eqs. (60/61). The wall region is too small (2 x 1073 GeV~), and is only in the 
figure For the chosen parameters, the ratio of the lepton density to the entropy density 
at is given by 
L3/s ~ 1074, (75) 


which results, after the inclusion of the dilution factor, Eq. (12) in a final baryon asym- 
metry after the electroweak phase transition of approximately 


ie eae (76) 
S 
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Figure 4: Particle number densities normalized to entropy as a function of spatial position 
z for a bubble whose wall is at z ~ 0 and parameters as described in the text. From top 


to bottom, the curves represent the following densities: h, Q3, br, L3, and tp. 
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exactly as observed. Of course, the particular value is highly dependent on our choice of 
parameters, but the ability of the Topflavor model to produce this value is not; the fact 
that the order of magnitude comes out correctly is indicitive of the fact that for natural 
values of parameters, this model can produce an appropriate baryon number. 

From our numerical integration of the differential equations, we can relax our equilib- 
rium assumptions and examine the effect of finite values for the rates, particularly I1. In 
Figure [5] we present the final lepton density as a function of T; (artificially assuming that 
the C'P-violating source, bubble parameters, and critical temperature are unchanged). 
For small T4, we see that the final lepton asymmetry is linearly proportional to T1. In 
this regime, the generated Higgs/top number densities are a constant background source, 
only a small fraction of which is converted into leptons by Ty. For large T4 instead, 
the equilibrium condition is reached and the dependence on I; saturates. We see from 
Figure [5]that our chosen parameters are just at the turn-on of this saturation region. 


4 Discussion 


We have explored the possibility that the baryon asymmetry of the Universe might have 
come from new gauge dynamics at the TeV scale. The primary impediment to realizing 
this idea is the washout by ordinary EW sphalerons, which we have avoided by having new 
physics which coupled asymmetrically to the three (lepton) generationd9, In particular, 
we showed that the Topflavor model, with SU(2); x SU(2)2 breaking at the TeV scale, 
can easily produce an appropriate baryon asymmetry for natural values of its parameters. 

It is interesting to compare the results obtained in Eq. (75) with the much smaller 
ones that are obtained by similar methods in the MSSM [19]. The difference stems from 
three main different factors: First, the magnitude of I; is much larger than the EW 
sphaleron rate Pws ~ 1074 GeV. As shown in Fig. [5] this induces an enhancement of a 
few orders of magnitude in the final result. Second, the value of A ~ 0.2 is two orders 
of magnitude larger than the typical values of AG ~ O (107°) obtained in the MSSM for 
moderate values of the CP-odd Higgs mass. Third, in the MSSM there is a suppression 
proceeding from the effective cancellation of chiral charges discussed in Refs. [54] [87] that 
is not present in the model analyzed in this article. 

We expect that the LHC can extensively test this idea. First, by discovering the 
SM Higgs boson and searching for the existence of other elements strongly coupled to it, 
the LHC can help construct the picture of the ordinary EW phase transition as first or 
second order. This would, for example, confirm or rule out the MSSM scenario of EW 
baryogenesis. Second, by discovering new elements, in the case of Topflavor, the heavy W’ 
and Z’ bosons with masses of about 2 TeV and preferred decays into third family fermions, 
should be easily discovered at the LHC [55]. Finally, by determining the properties of 
the new elements, one may determine whether or not they are viable as a mechanism of 


6A parallel idea, which we chose not to explore in detail, is to have new physics coupled asymmetrically 
between quarks and leptons, as in [53]. Anomaly cancellation is very challenging in this framework, and 


the resulting models are therefore somewhat contrived. 
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Figure 5: Lepton number density normalized to entropy as a function of Ty. The red 
dashed line is the bound on I; from non-observation of proton decay mediated by instan- 
tons, whereas the blue dotted line is the bound inferred by requiring that the broken phase 
minimum remain the true vacuum at T = 0. The star indicates the sample parameter 


point considered in the text. 
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baryogenesis. In the case of Topflavor, the need for a first order phase transition in the 
breakdown of SU(2), x SU(2)z is connected to both a small quartic for the Higgs o and a 
reasonably strong gauge coupling gı. The first of these has the dramatic consequence that 
the mass of ø is well below the symmetry-breaking scale, m? ~ \u?. For the parameters 
we explore in detail, m, ~ 600 GeV, making o the lightest of the new particles. It couples 
strongly to the Higgs sector, and to the new gauge bosons. 

An extended gauge sector at the TeV scale has motivations from many mysteries in 
particle physics. In addition, it allows a much less constrained structure to explore the 
idea of baryogenesis through gauge symmetry breaking. The SM picture arising from the 
EW breakdown is economical, but is under assault from the null LEP Higgs searches, 
favoring a second order phase transition, and from the lack of sufficient C P violation to 
produce enough asymmetry. Alternatives based on the ordinary EW transition such as in 
the MSSM can be viable, but remain constrained by the Higgs mass, and from the fact 
that the new sources of C'P violation remain tightly connected to low energy C'P-violating 
observables. Extended gauge symmetries can naturally have first order transitions, and 
CP violation is further removed from low energy observables. The primary obstacle is 
the fact that the EW sphalerons will try to erase any generated baryon asymmetry with 
B — L = 0, but this is overcome if the SM matter does not couple universally to the new 
dynamics. Ultimately, the LHC will help in resolving this question. 
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